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We study the evolution of the single-hole spectral function when the ground state of the 
anisotropic-triangular-lattice antiferromagnet changes from the incommensurate magnetically- 
ordered phase to the spin-liquid state. In order to describe both of the ground states on equal 
footing, we use the large-A'' approach where the transition between these two phases can be ob- 
tained by controlling the quantum fluctuations via an 'effective' spin magnitude. Adding a hole 
into these ground states is described by a t-J type model in the slave-fermion representation. Im- 
plications of our results to possible future ARPES experiments on insulating frustrated magnets, 
especially CS2CUCI4, are discussed. 
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I. INTRODUCTION 

The physics of geometrically frustrated quantum mag- 
netic systems has been an arena for discovery of new and 
interesting phenomena both on the theoretical and ex- 
perimental fronts. Much of the effort has been sparked 
by Anderson's proposal of the exotic resonating valence 
bond (RVB) spin-liquid phasei*^ and the belief that frus- 
trated systems are likely to possess RVB spin-liquid-like 
ground states in two dimensions. Theoretical efforts in 
the past have explored a variety of two-dimensional frus- 
trated systems^''*'^. The discovery of superconductivity^ 
in Nao.35CoO2-l.3H2O has also raised the hope that the 
RVB could be the cause of superconductivity in certain 
frustrated magnetic systems. 

Recent theoretical works on the insulating anisotropic- 
triangular-lattice antiferromagneliiSiSiiSiii have been 
particularly motivated by the experimental discovery 
of the fractionalized spin-liquid phase in Cs2CuCl4i24i^, 
where spin-1/2 moments of Cu antiferromagnetically in- 
teract with each other in the layers of anisotropic trian- 
gular lattices. The experiment has shown that, with a 
sufficiently high magnetic field within its ordering plane, 
the compound looses its long-ranged cycloid order and 
enters the spin- liquid phaseiSiiS. fn this phase, the dy- 
namical spin structure factor shows a broad profile char- 
acteristic of fractionalization of spin-1 magnons into a 
pair of spin-1/2 spinonsi^ii^. 

A useful description of this transition can be obtained 
in the large- TV approach^d'^d'^d^d^d^d^i^^i^^ to the anti- 
ferromagnetic Heisenberg model on the anisotropic tri- 
angular lattice. Here the spin SU(2) symmetry is gen- 
eralized to Sp(iV) and the large-A'' limit of the model 
is investigated. In this approach, the quantum fluctua- 
tions are controlled by the parameter k that is the Sp(Af) 
analog of the spin magnitude. The resulting phase dia- 
gram contains both the incommensurate long-ranged or- 
der (LRO) phase^^ (this corresponds to the cycloid order) 
and the spin-liquid stateS^ in the large and small k limits, 
respectivelySii^. Thus one may expect that 1/k plays a 



role similar to the in-plane magnetic field in CS2CUCI4. 
Indeed the spin structure factor of the spin-liquid state 
computed in the large- theorjiSiiS looks very similar to 
that observed in the neutron scattering experiments on 
Cs2CuCl4i2Ja. 

To this date, experimental evidence of the RVB spin- 
liquid phase in the anisotropic triangular lattice has 
been obtained only in neutron scattering experiments on 
Cs2CuCl4''^^^»i^. Thus it is desirable to have independent 
experimental confirmations of these results. One such ex- 
periment may be the angle-resolved photoemission spec- 
troscopy (ARPES) that measures the electron spectral 
function. In particular, the result of ARPES on an insu- 
lating compound at half-filling corresponds to measuring 
the single-hole spectral function in a given ground state 
of the insulator—. 

In this paper, we investigate the evolution of the single- 
hole spectral function across the transition between the 
incommensurate LRO phase and the spin-liquid phase. 
Previously the single-hole spectral function in the LRO 
phase was investigated by taking into account the in- 
teraction of the hole with the spin wave32S*2SiSL2&. The 
well-defined low-energy quasiparticle peak was found in 
these studies. This is due to the fact that small low- 
energy density of states of the magnons leads to only a 
weak perturbation to the quasiparticle coherence at low 
energy. While this calculation is enough to describe the 
spectral function deep inside the ordered phase, an addi- 
tional ingredient may have to be considered close to the 
transition toward the nearby spin-liquid state where the 
spinous are deconfined— . The confinement of the spinous 
in the LRO phase is a low energy phenomenon and it is 
expected that the confinement energy scale will be very 
small near the transition. Therefore, it is expected that 
there will be spinon excitations at intermediate to high 
energies near the transition and they may affect the co- 
herence of the hole motion. For example, the neutron 
scattering experiments on the magnetically ordered phase 
of CS2CUCI4 found the spinon continuum as well as the 
magnon peak, implying that this ordered phase is very 
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close to the transitioniSii^. Indeed the dominant inter- 
action near the transition may be the spinon-hole inter- 
action since the magnons are relatively suppressed close 
to the transition. 

We use the results of the large- iV theory to describe 
the transition between two ground states at half-filling. 
In this formulation, as will be explained later, the 
spinous appear at the mean field level even in the LRO 
p]^j^gg8ji4ji5ji6ji7^ Even though this is an artifact of the 
mean field theory and the low energy spinous are con- 
fined due to the fluctuation effects beyond the mean field 
thcoryMi^^t^L^, we can still use these spinous at inter- 
mediate to high energies to mimic the situations near the 
transition as described above. The motion of the single- 
hole is incorporated via a t-J type model in the slave 
fermion representation^^. Our computation takes into 
account the multiple interactions between the motion of 
the hole and the spin background, but does not include 
the vertex corrections. In principle, the latter approxi- 
mation could lead to an underestimation of the coherence 
in the motion of the hole. Comparison to the exact di- 
agonalization studySS, however, suggests that the latter 
approximation does not change the overall structure of 
the spectral function in a significant fashion. 

Our results indicate that the single-hole spectral func- 
tion still shows a well-defined low-energy quasiparticle 
peak in the incommensurate LRO phase with a visible 
dispersion in momentum space. This is due to the fact 
that the density of states of spinons becomes smaller as 
the energy scale is decreased and, as a result, does not 
strongly disturb the low energy motion of the hole. As 
the transition to the spin-liquid state is approached, the 
spectral weight of the quasiparticle peak is decreased and 
the incoherent continuum at intermediate to high ener- 
gies grows. Eventually the quasiparticle peak disappears 
at the transition. In the spin-liquid state, the single-hole 
spectral function only shows a broad incoherent contin- 
uum with little dispersion, which can be interpreted as a 
consequence of spin-charge separation. 

The rest of the paper is organized as follows. In sec- 
tion II, we describe the formalism used to compute the 
single-hole spectral function and explain the approxima- 
tion schemes used in this work. In section III, the numer- 
ical results on the single-hole spectral function in various 
phases are discussed. Finally, we summarize our results 
and conclude in section IV. 



II. FORMULATION OF THE PROBLEM 

At half-filling, the quantum ground states of the Mott 
insulator on the anisotropic triangular lattice are de- 
scribed by the Heiseuberg model: 

ifj = Ji ^ S, • + J2 ^ S, • S,, (1) 

{i,3) iihj)) 

where the first and second terms represent the nearest 
and the next-nearest neighbour Heisenberg couplings, re- 
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FIG. 1: Anisotropic triangular lattice with two types of bond. 
The lattice is modelled by a square lattice with an additional 
interaction along one of the diagonals. The nearest neighbour 
coupling is associated with subscript 1 and the next-nearest 
neighbour coupling is associated with subscript 2. 



spectively (see Fig.l). The motion of the single hole in a 
given ground state of the Mott insulator can be studied 
by adding the kinetic energy terms of the electrons; 

iit = -ti + /i.e.) -t^Y. (^l^.'- + '^•^■)' (2) 

where cj^^. is the creation operator of the electrons with 
spin CT =t, I (repeated indices are summed) and it is as- 
sumed that there is no doubly occupied site in the Hilbert 
space to take into account the strong correlation inherent 
in the Mott insulator. 



A. Ground States of the Mott Insulator via the 
Large-A Sp(A) Heisenberg Model 

A useful formulation of the Heisenberg model in EqQ] 
can be obtained from a bosonic representation of the 
SU(2) spin via 5*^ = \h\^T''^^,h„, , where r'' (a = x,y,z) 
are the Pauli matrices and hi„ is the canonical bosonic 
operator. In addition, the constraint bl^bicr = rib = 2S 
must be imposed at each site, where S is the magnitude 
of the spin. In order to study this model on frustrated 
lattices in a controlled fashion and to obtain results that 
are non-perturbative in the coupling constants, we con- 
sider the Sp(A^) generalizationi^ii^iiSiil of the physical 
spin SU(2) ~ Sp(l). This can be achieved by introducing 
N flavours of bosons at each site, bia, where a = (n, a) 
and n = 1, is the flavour index. Here again the bo- 
son number at each site, rif, = bj^bia, must be fixed. The 
resulting Sp(A^) Heisenberg model is given 

HspiN) - - E ^i'^c.0blbl){J^sb^,b,s), (3) 

i>j 

where = Ji and J2 on the nearest and next-nearest 
neighbour links, respectively. Here J'ap = —Jpa is a 
generalization of the SU(2) antisymmetric tensor Co-o-' = 
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FIG. 2: The zero temperature phase diagrams' of the Sp(A^) 
Heisenberg antiferromagnet on the anisotropic triangular lat- 
tice at large A^. 1/k represents the strength of quantum fluc- 
tuations. Single hole spectral functions are evaluated at each 
of the ground states denoted by the circles along paths (1) 
and (2). See text for details. 



—ta'a and is given by a 2N x 2N matrix that has N 
copies of e along its center block diagonal. The matrix is 
zero elsewhere. 

We consider the large- iV limit of this model at a fixed 
boson density per flavour, k — Ub/N . Notice that k — 2S 
when = 1. Thus k is the Sp(iV) analog of 2S and 
controls the strength of quantum fluctuations. The large- 
N mean-fleld theory of this model has been analyzed^^ 
and the phase diagram at large- has been obtained as 
a function of k and J2/{Ji + J2) (see Fig. 2)^. For large 
values of k, quantum fluctuations are small so that the 
LRO states arise by breaking the global Sp(A^) symmetry. 
In the small k limit, strong quantum fluctuations lead to 
paramagnetic states with only short-ranged order (SRO). 
In the phase diagram of Fig. 2, the ordering wavevector of 
the LRO states (and the corresponding SRO phases) are 
labeled by the ordering wavevector q = (91, 92) in two 
dimensions. It has been known that (tt, tt) SRO phase 
develops the valence-bond-solid order upon inclusion of 
singular fluctuations beyond the mean-field theory and 
the Berry phase effect^. In this paper, we will consider 
the transition between the incommensurate {q, q) LRO 
and the (g, q) SRO phases that respectively correspond 
to the cycloid-ordered and the translationally-invariant 
spin-liquid phases discovered in the experiment involving 
Cs2CuCl^i2*i^. 

Now, for future reference, some explanation of the 
large-TV mean-field theory is in order. The presence of 
the LRO can be described by the condensation of bosons, 
(6"'^) = -\/iV<5"a;,f, where the c-number field xf quantifies 
the condensate. Thus, in the LRO phase, it is convenient 



to parametrize the bosonic fields as 



(4) 



where m = 2,...,N. On the other hand, none of the 
bosons is condensed in the SRO phase. After decoupling 
the quartic boson interaction in Eq|3| by the Hubbard- 
Stratonovic fields Qij, the mean- field Hamiltonian can 
be obtained asiii^JMiiia 



Hmf = 



i>j 



N 



J ij 



N 
m=2 



+ ^A, |iV|<|^ + 



N 



rib 



(5) 



where Qij = Tq{Japb\j)jp) at the saddle of the corre- 
sponding action. Here the last term of Eq|Sltakes care of 
the constraint on the boson number at each site and Ai 
is the Lagrange multiplier. 

At the mean- field level, Qij = Qi and Q2 on the near- 
est and the next- nearest neighbour links, and A is a site- 
independent constant. In the LRO phase, the boson con- 
densate fraction has the form 



xe 



iQri 



-IX e 



(6) 



where q = 2Q is the magnetic ordering wavevector. In 
the mean-field theory, Qi, Q2^ x, A and Q have to be de- 
termined self-consistcntly for a given k — Ub/N . 

The part of the Hamiltonian for the uncondensed 
bosons or spinous, foimo- (w = 2, ...,7V), can be diago- 
nalized using the Bogoliubov transformation: 

N 

H'mF = mi ("km "km + /^km^m) , (7) 



k m=2 



where 



7k = i [JiQi(sinfea; -I- sinfcj/) 



J2Q2sin(A:3, -I- ky)] 



(8) 



The Bogoliubov quasiparticle operators, akm and /3km, 
are given by 



akm 

^-km 



Uk l^k 
Wk 



0km T 
^-kmi 



where 



Uk 



1 

72 



A + I, 

Ck 



Wk 



V2 



A 

1. 

Ek 



(9) 



(10) 



Notice that, in the SRO phase, there are no condensed 
bosons and all the spinous are described by Eq[7| 
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It has been known that the presence of unconfined 
spinons is an artifact of the large-A^ mean-field the- 
ory and these spinons are, in principle, confined at 
low energies due to fluctuations beyond the mean-field 
theorjiMiiSiiSJ'^ Thus the spinons cannot be seen deep in- 
side the LRO phases. The confinement scale below which 
the confinement occurs, however, can be very small when 
the transition to the spin-liquid phase is approached. 
Thus we expect that there exist 'almost'-deconfined (or 
loosely-confined) spinons at intermediate to high energy 
scales near the transition. At the transition point, this 
energy scale should be zero because the spinons are de- 
confined in the spin-liquid phase. In fact, the experi- 
ment on CS2CUCI4 reveals the existence of both of the 
magnon peak and the spinon continuum in the cycloid- 
ordered phase, indicating that this system is very close 
to the transition toward the spin-liquid phaseiSii^. While 
the confinement phenomenon cannot be captured by the 
mean-field theory, the existence of the deconfined spinons 
at the mean-field level can be used to phenomenologically 
describe the existence of 'almost'-deconfined spinons near 
the transition. From now on, we will adopt this picture 
and examine the effect of the 'almost'-deconfined spinons 
on the motion of the holes. 



B. Motion of Single Hole and the Spectral Function 

In order to describe the motion of a single hole in the 
large- limit, we use the slave- fermion representation of 
the electron operator with the flavour index: 

"^l,Ti<T ^ fi ^l,ncn (11) 

where fi represents the fermionic charge degree of free- 
dom and bi^na the bosonic spin degree of freedom of the 
electrons. 

The single-hole spectral function can be computed 
from the single-hole Green's function: 

Gh,a{r,T) = ^^Gh,na{r,T) 

n 

= -ljT.('^rcUr,r)cnAO,0)), (12) 

n 

where T-r is the time-ordering operator of the imaginary 
time. Notice that the physical hole Green's-function is 
defined as the average of the Green's functions in differ- 
ent flavour channels. Since no flavour degree of freedom 
is special on physical grounds, this is a legitimate defi- 
nition of the hole Green's function that turns out to be 
particularly useful for our purposes. 

After analytic continuation, the spectral function can 
be obtained from 

AhA^,Lo) = --lm[GhM^,u + tS)l (13) 

where G/i,(j(k, cj) is the single- hole Green's function in 
Fourier space. 



At the mean-field level, the single-hole Green's func- 
tion can be written as 

GnAr. ^) « E ^/('•' ^)Gb,nAr, r), (14) 

n 

where 

G/(r,T) = -(T,/(r,r)/t(0,0)) 
G6,„^(r,r) = -(T,6L(r,r)&„.(0,0)). (15) 

In Fourier space, we get 
G'ft,o-(k,Lj) 

==-^EX] f di^ Gb^n„{k~ q,uj ~ i^)Gf{<i,i^). (16) 

Thus the single-hole Green's function is given by the con- 
volution of the fermionic and the bosonic Green's func- 
tions. 

In this work, we will assume that the injection of the 
singic-holc docs not affect the underlying ground state at 
half-filling. In this case, it may be sufficient to use the 
bare bosonic Green's-function, G°„j^, in the absence of 
the hole^^. On the other hand, the fermionic degree of 
freedom of the hole will be very much affected by the pres- 
ence of the excitations in the underlying ground state. 
We will compute the renormalized fermionic Green's 
function using the self-consistent Dyson equation2L2&, 
where the non-crossing diagrams are summed. Possible 
effect of the vertex corrections in the fermionic Green's 
function and EatTHwill be discussed later. 

1. Computation of the Fermion Green's Function 

The proper large- A'^ generalization of the kinetic energy 
terms of the holes in Eq|2| is given by 

^* = -JfT.(f4n.hn.f]+h.C.) 

- I E ifM,n.hn<rf^, + h.C.) (17) 

Fourier transformation of these terms leads to 

- E Tk-k'/k^l, „,fok'-q,„./Lq (18) 

kk'q 

where 

Fk = 2 [ti (cos kx -\- cos ky) + t2 cos{kx + ky)] . (19) 

These hopping terms give rise to the interaction between 
the hole and the spin background. 

In the LRO phase, the existence of the condensed 
bosons, |(6k,ncr)P = A^^r|a;p^k,Q, leads to 

= 2^rk-Q|a:|2/t/k 

k 



FIG. 3: The self-energy correction of Eg 1211 The wavy hne 
represents the fermionic propagator and the sohd hne repre- 
sents the bosonic propagator. 
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-Q,m<T flh-ci (20) 



kk'q m=2 



The second term corresponds to the interaction between 
the uncondensed spinons and the fermion and gives rise 
to the fermion self-energy corrections. 

Let us first consider the perturbative self-energy cor- 
rections in ti/N and t2/N. The lowest-order correction 
is given by (see Fig. 3) 



Si(k) = 2^rk_k'kk' 



(21) 



where {N — 1)/A^ « 1 is used. The next order correction 
involves the diagram in Fig. 4 and can be computed as 

E2 (k, c^) = i- ^ /kk'qG/ {k~q,u;- Ek' - eq-k' ) , (22) 



k',q 



where 



/kk'q = iFk-k'^k'^^k'-q + Fk+k'-qUk'-q^'k' 



(23) 



and (TV - 1)/N^ « 1/A^ is used. Notice that, in 
the perturbative calculation, the bare fermionic Green's- 
function, Gj, should be used for G/ in Eq[221 

In order to take into account the multiple interactions, 
however, we need to go beyond the perturbative evalu- 
ations of the self-energy contributions. In our work, we 
will determine the renormalized fermion Green's func- 
tion, G/, by summing up diagrams in the non-crossing 
approximation. This amounts to using the renormal- 
ized Green's function in Eg 1221 and solving the following 
Dyson equation for the fermion Green's function. 



Gj\k,Lu) =uj + 2|:Eprk-Q - I]i(k) - I]2(k, 



(24) 



In the numerical solution of this equation, we set the 
coefficient of the self-energy correction in Eq|221to be one 
for convenience {i.e. 1/N = 1). Notice that 2|a;pFk-Q — 
Si(k) corresponds to the dynamically generated 'band' 
dispersion for the fermions at the lowest order. In the 
case of the SRO phase, the same Dyson equation can be 
used with a; = in Eg OH 



FIG. 4: The self-energy correction of Eg 1221 The thick wavy 
line is the exact fermionic propagator. Combined with the 
contribution from Fig|3| this gives rise to the Dyson equation 
for the fermionic propagator in Eg 1241 



2. Evaluation of the Single-Hole Spectral Function 

In the LRO phase, the Green's function for the con- 
densed bosons in the Fourier space is given by 

Gb,iT(k,c^) - -^N\x\^Sk.^QSiLo) 
Gb,u{k,Lo) - -N\x\''SumS{uj). (25) 

On the other hand, the Green's function for the uncon- 
densed bosons are obtained as (for both a =t, j) 



Gb^mo- (k, + iS) = 



\vk\ 



(oj - Ek + iS) (oj + Ck + iS) 



(26) 



Then the single-hole Green's function in Eq^] can be 
written as 

N - 1 

Gh,a = Gh,la H — Gh,2(T ~ Gh^la + Gk,2a- (27) 

The single-hole spectral function, A^^a — — (l/7r)ImG?i,CT 
is now given by 

A,,,^(k, uj) = Ah,ia(M, ^) + A,2^(k, uj). (28) 
The first contribution from the condensed bosons is 

A/i,it(k,w) = --\x\'^lmGf{k + Cl,uj + id), 

TT 

Ahu(k,uj) = --\x\^lmGf{k-Cl,LU + iS). (29) 

TT 

The second term is due to the interaction between the 
fermion and the uncondensed bosons and is given by 

Ah^2aik,Uj) 

= - ^ / dyAf (q, y)Ab^2a (q - k, y - cj) (tj > 0) 
q -'o 

^Y. dyAf{ci,y)AiaM-Kv~^) (^ < 0), (30) 

q 

where Af = — (l/7r)ImG/ and A\j,2a = — (l/7i')IinGfc^2(T- 

On the other hand, in the spin-liquid state (or the SRO 
phase), the first term in Eg 1281 vanishes because there is 
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FIG. 5: Eight momenta in the principal Brillouin zone at 
which the spectral functions are evaluated. 



no boson condensate and the spectral function is entirely 
given by the second term. 

We will first numerically solve the Dyson equation for 
the fcrmion Green's function in both the LRO and the 
spin- liquid phases. Then the fcrmion spectral function 
can be convoluted with the boson spectral functions in 
the two phases to get the expression for the single-hole 
spectral function. The evolution of the resulting single- 
hole spectral function across the transition between the 
LRO and the spin-liquid phases is discussed in the next 
section. 



III. NUMERICAL RESULTS 

At first, the fermion Green's function is obtained by 
solving numerically the self-consistent Dyson equation on 
a discrete mesh of k and uj points. In order to do so, 
we have solved the saddle point equation of the large- TV 
mean- field theory for each set of k, Ji, J2, and use the 
mean field results as the input to the Dyson equation 
of the fermion Green's function. The full hole Green's- 
function is then computed by using Eas l28l29l30l The 
calculations are done with a finite lattice of size 8^. We 
have confirmed that the overall structure of the spectral 
function does not change when the size of the system is 
increased to 16^. 

The evolution of the spectral function near the tran- 
sition between the incommensurate LRO and the spin- 
liquid phases, is examined at 8 different points in the 
phase diagram by following two different paths (see 
Fig|2)- On both paths, the transition between two states 
is achieved by changing k for a given J2/(Ji -I- J2)- For 
each ground state, the spectral function is calculated at 
eight different momenta in the principal Brillouin zone. 
These eight momenta are shown in Fig|S| When the spec- 
tral functions are plotted for these eight momenta, they 
will be arranged from bottom up in the order indicated 
by the bracketed numbers as shown in FigEl All other 
similar plots will follow the same convention. 

As explained earlier, quantum fluctuations are stronger 
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FIG. 6: Evolution of the spectral functions along path 1 in 
Fig. H Ji = J2 = 1.0 and U/Ji = 12/ J2 = 2.5. 



for smaller values of k. Thus, the transition from the 
incommensurate LRO phase to the spin-liquid state can 
be obtained by decreasing k for a given value of J2/{Ji + 
J2). We consider paths (1) and (2) corresponding to 
J\l J2 — 1 and J1/J2 = 0.75, respectively. The results 
for the four different values of k are shown in Fig El for 
path (1) and in Fig[7|for path (2). 

The results of paths (1) and (2) are very similar. At 
large k, there exists a sharp low-energy peak in the in- 
commensurate LRO phase suggesting the existence of a 
coherent hole excitation or a quasiparticle. Notice that 
the incoherent continuum at higher energies are due to 
the existence of the 'unconfined' spinous. In the spin- 
liquid phase, the quasiparticle peak completely disap- 
pears and only the incoherent broad continuum can be 
seen in the spectral function. 

Notice that the spectral weight of the quasiparticle 
peak in the LRO phase decreases as the transition to the 
spin-liquid phase is approached and instead the spectral 
weight of the incoherent continuum grows. This is due 
to the fact that the spectral weight of the quasiparticle 
peak is proportional to the strength of the boson conden- 
sate or the magnetic order parameter that decreases as 
the transition point is approached. This can be observed 
in both FigEland Fig[7| but is illustrated more clearly in 
the latter figure. In Fig[7| the strong quasiparticle peak 
seen when k = 1 diminishes as we approach the transi- 
tion. The remnant of the peak is indicated by the arrows 
in the figure. Upon careful inspection, a slightly broad 
bump beside the peak consistently persists in all of the 
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FIG. 7: Evolution of the spectral functions along path 2 in 
Fig. H Ji = 0.75, J2 = 1.0 and fi/Ji = f2/J2 = 2.5. 



plots indicating that it is a contribution from the uncon- 
fincd spinons. One can observe in both figures that as 
the spectral weight of the quasiparticle peak shrinks, the 
remaining spectral weight gradually moves to the 'uncon- 
fined' spinons at higher energies. 

The incoherent continuum in the SRO phase arises be- 
cause the injected hole decays into its charge and spin 
degrees of freedom and they do not propagate together 
in a coherent fashion. The origin of the broad continuum 
in the LRO phase is basically the same; the 'unconfined' 
spinons make the motion of the hole incoherent. The 
reason why the low-energy quasiparticle peak survives in 
the LRO phase is that the density of states of the spinons 
becomes very small in the low energy limit so that the 
spin degree of freedom of the injected hole at low energy 
cannot efficiently decay into the spinon continuum. 

We have also confirmed that varying the hopping am- 
plitude, t, does not change the qualitative behaviour of 
the spectral function. The only change is that the band- 
width of the quasiparticle is increased as the hopping 
amplitude becomes bigger. 



IV. SUMMARY AND CONCLUSION 

We investigated the single-hole spectral function in 
the anisotropic-triangular-lattice Heisenberg antiferro- 
magnet. In particular, we studied the evolution of the 
single-hole spectral function across the transition be- 
tween the incommensurate LRO phase and the spin- 



liquid state. This phase transition can be described by 
the large- iV Sp(A^) mean- field theorjiiliSiiSiiLia, where 
the Sp(A^) analog of the spin magnitude, k, can be ad- 
justed to control the strength of quantum fiuctuations. 
The motion of the single hole injected into the insulator 
is studied via the t-J type model in the slave fermion rep- 
resentation, where the electron is represented as the com- 
posite of a charge-carrying fermion and a spin-carrying 
boson. Taking into account the multiple interactions be- 
tween the injected hole and the excitations of the un- 
derlying insulator in the non-crossing approximation, we 
evaluated the single-hole spectral function in both of the 
LRO and the spin-liquid phases. 

It is found that, in the LRO phase, the quasiparticle 
peak exists at low energy and the broad incoherent back- 
ground also arises at intermediate to high energies. This 
incoherent continuum arises due to the existence of 'un- 
confined' spinons in the large-TV mean-field theory. In 
principle, these spinons should be confined below a 'con- 
finement' energy scale due to a fluctuation effect beyond 
the mean-field theory^. Near the transition, however, the 
confinement energy scale becomes very small (in fact it is 
zero at the transition) so that the spinons at intermediate 
to high energies are 'almost' unconfined. In this respect, 
our calculation is particularly suited to describe the spec- 
tral function near the transition. As the transition to the 
spin-liquid state is approached, the spectral weight of the 
quasiparticle peak goes to zero and the incoherent con- 
tinuum grows. Finally, in the spin-liquid phase, only the 
incoherent continuum survives in the spectral function. 
This is due to the fact that the injected hole can de- 
cay into the charge-carrying fcrmionic excitation and the 
spin-carrying spinon and the motions of these excitations 
are not coherent together. 

In our studies, we neglected various vertex corrections, 
which may have lead to an underestimation of the co- 
herence, especially in the spin-liquid phase. However, 
the recent exact diagonalization study of the spin-liquid 
state on the Kagome lattice showed spectral functions 
that were indeed completely incoherent, leading to the 
conclusion that the vertex correction would not change 
the qualitative behaviours obtained in our study. 

Our work is relevant to possible ARPES experiments 
on Mott insulators on the anisotropic triangular lattice 
since the ARPES on the insulator can measure the single- 
hole spectral functionSi. In particular, notice that the 
neutron scattering experiment on CS2CUCI4 reveals the 
existence of the spinon continuum at intermediate to high 
energies in the LRO phase near the transition as well as 
in the spin- liquid phas o^^i^^ . The parameter 1/k in the 
large- TV theory plays the same role as the in-plane mag- 
netic field in the experiment on CS2CUCI4, where the 
transition between the two phases is achieved by chang- 
ing the in-plane field or temperature. Thus we expect 
that, in light of our work, the ARPES experiments on 
such systems will provide independent confirmation of 
the spinon continuum discovered in the neutron scatter- 
ing experiments on CS2CUCI4 or in related systems. 
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